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Abstract 


We present recent improvements of the modeling of the disruption of strength dominated 
bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements 
include an updated strength model and a friction model, which are successfully tested by 
a comparison with laboratory experiments. In the modeling of catastrophic disruptions of 
asteroids, a comparison between old and new strength models shows no signihcant deviation 
in the case of targets which are initially non-porous, fully intact and have a homogeneous 


structure (such as the targets used in the study by Benz and Asphaug 1999). However, 
for many cases (e.g. initially partly or fully damaged targets, rubble-pile structures, etc.) 
we hnd that it is crucial that friction is taken into account and the material has a pressure 
dependent shear strength. Our investigations of the catastrophic disruption threshold as 
a function of target properties and target sizes up to a few 100 km show that a fully damaged 
target modeled without friction has a which is signihcantly (5-10 times) smaller than 
in the case where friction is included. When the effect of the energy dissipation due to 
compaction (pore crushing) is taken into account as well, the targets become even stronger 
{Q*j^ is increased by a factor of 2-3). On the other hand, cohesion is found to have an 
negligible effect at large scales and is only important at scales < 1km. 

Our results show the relative effects of strength, friction and porosity on the outcome 
of collisions among small (< 1000 km) bodies. These results will be used in a future study 
to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt and Stewart 


2012 ) 


Keywords: Asteroids, collisions, Collisional physics 


1. Introduction 

Collisions play a fundamental role in the formation and evolution of the planets and the 
small body populations of the Solar System. Models of the evolution of such populations 
(e.g. the Asteroid Belt) compute the time dependent size and velocity distributions of the 
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objects as a result of both collisional and dynamical processes. A scaling parameter often 
used in such numerical models is the critical specific impact energy which results in the 
escape of half of the target’s mass in a collision. The parameter is called the catastrophic 
impact energy threshold (also called the dispersion threshold). The specific impact energy 
is often defined as Q = Q.brripVp/M t, where rup, Vp and Mt are the mass and speed of the 
projectile and the mass of the target, respectively. The catastrophic disruption threshold 
is then given by the specific impact energy leading to a largest (reaccumulated) fragment 
containing 50% of the original targets mass. In recent studies (e.g. Stewart and Leinhardt 


2009 Leinhardt and Stewart 2012), a more general definition of the specihc impact energy 


was proposed which also takes the mass of the impactor into account: 


Qr — 


0.5mpVp + 0.5MtV^ 0.5/xld^ 


Mi 


tot 


Mi 


( 1 ) 


tot 


where Mtot = rnp + Mt and /ij = rUpMT/Mtot and Vi is the relative velocity. The corre¬ 
sponding radius Rci is defined as the spherical radius of the combined projectile and target 
masses at a density of 1 g cm“^. According to this new definition, the catastrophic disruption 
threshold is then called 

Values of Q*^, (or Q*jit>) have been estimated using both laboratory and numerical hy¬ 
drocode experiments (see e.g. Holsapple et ah 2002 Asphaug et ah 2002). For the small 


body populations, the hrst suite of numerical calculations aimed at characterizing the catas¬ 
trophic disruption threshold in both the strength regime and the gravity regime was per- 


formed by Benz and Asphaug 

(1999 

, who used a smoothed particle hydrodynamics (SPH) 

code ( 

Benz and Asphaug 

1994 

1995 

) to simulate the breakup of basalt and icy bodies from 


centimeters-scale to hundreds kilometers in diameter. More recently, [Leinhardt and Stewart 
(2009) computed Q}) curves using the hydro code CTH (McGlaun, 1990) to compute the 


fragmentation phase and the N-body code pkdgrav (Richardson et ah, 2000) to compute the 


subsequent gravitational evolution of the fragments. In this study, the dependency of Q*jy 
on the strength of the target was investigated. In a recent study by Jutzi et ah (2010), the 


effect of target porosity on was investigated using an extended version of the SPH code 
(Jutzi et ah, 2008). In this study, the size and velocity distribution of the fragments was 


computed as well, using the pkdgrav code. Benavidez et ah (2012) performed a study of a 


large number of collisions among Rt = 50 km rubble pile bodies using the original SPH code 


by Benz and Asphaug (1994, 1995). 


A numerical tool which is very suitable and has been often used to study disruptive 
collisions among rocky bodies in general, and was used in many of the above-mentioned 
asteroid disruption studies, is based on the SPH method. Over the last decades, the basic 
method has been extended by implementing additional physics (e.g. Benz and Asphaug 


1994, 1995 Jutzi et ah 2008, 2009) with the goal to realistically model rocky bodies with 


various internal structures. In addition to improved constitutive models, models which 
mimic the complex macroscopic structure of rubble pile-like bodies have been used as well 
(e.g. Asphaug et ah 1998 Benavidez et al. 2012). Although the SPH models used in 
planetary sciences have been signihcantly improved over the last decades, they were still 
lacking aspects that can be important in some impact regimes. For example, previous SPH 
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strength models used the so-called von Mises yield criterion, which does not describe well the 
behavior of rocky materials which are known to have a pressure dependent shear strength. 
Furthermore, fully damaged material was treated as a strengthless ’fluid’ in previous models. 
Finally, while self-gravity is included in the versions of the SPH codes used to model giant 
collisions, it is often not implemented in the SPH code versions which also include the physics 
of solid bodies. 

In this paper, we present improvements of the SPH technique concerning the modeling 
of the disruption of strength dominated bodies. These improvements include a pressure 
dependent Drucker-Prager-like yield criterion, and a friction model for the damaged material. 
A few test cases are presented. Using the improved models, we then systematically study 
the effects various target properties (strength, porosity and friction) on the outcome of 
a disruptive collision (QJj). In this study, we use targets with a homogeneous internal 
structure (i.e., it is assumed that the voids or inhomogeneities are sufficiently small that 
their distribution can be assumed uniform and isotropic over the relevant scales). In section 
2, we present our numerical tool and the recent improvements and show two test cases. In 
section 3, the catastrophic disruption threshold Q*jj is investigated as a function of material 
properties. In section 4, the results are discussed and future work is indicated. 


2. Modeling 


2.1. Previous SPH models 


Benz and Asphau^ ( |1994| |1995[ ) extended the standard gas dynamics SPH approach to 

e.g. 


include an elastic-perfectly plastic material description 


see, 


1991) and a model of brittle failure based on the one of Grady and Kipp (1980). In the 


Libersky and Petschek 


fracture model, a state variable D (for damage) was introduced which expresses the reduction 
in strength under tensile loading and which varies between D = 0 and D = 1. Damage 
accumulates when the local tensile strain reached the activation threshold of a flaw. As 
stress limiter, the von Mises yield criterion was used. Finally the so-called Tillotson equation 


of state for basalt (Tillotson, 1962) was used to relate the pressure to the density and the 


internal energy. We refer the reader to the papers by 

Benz and Asphaug (1994, 

1995 

) for 

a detailed description of this method. This code was 

then used for instance by 

Benz and 


Asphaug (1999) to make a hrst complete characterization of for basalt and ice targets 
at different impact speeds. The same version (in terms of material models) of the SPH code 


was used by Benavidez et ah (2012) to study collisions among rubble pile asteroids. Jutzi 


et ah (2008) extended the Benz and Asphaug (1994, 1995) method by implementing a sub¬ 


resolution porosity model based on the P-alpha model (Herrmann, 1969 Carroll and Holt 


1972). In this implementation, a distention parameter a = ps/p is introduced where p is the 
bulk density of the porous material and pa is the density of the corresponding solid (matrix) 
material. The distention a is dehned as a function of pressure via the so-called crush-curve. 
It is used in the EOS to compute the pressure as a function of the matrix density. The 
porosity model was successfully tested by a comparison to laboratory experiments involving 
porous pumice (Jutzi et ah| 2009). The SPH code used by Jutzi et ah (2008) and in the later 
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(1994, 1995) 


studies is a parallelized version (Nyffeler, 2004) of the original code by Benz and Asphaug 


2.2. Recent improvements 

In the original implementation of the strength and fracture model by Benz and Asphaug 
(1994, 1995), the pressure independent von Mises yield criterion was used. However, it is 


known that the shear strength of rocks is pressure depended (i.e., it increases with increasing 
conhning pressure). A pressure dependent yield criterion often used to deal with rocky 
materials is the Drucker - Prager yield criterion 


>2 + Oifpli — kc — 0 (2) 

where R is the hrst invariant of the stress tensor, J 2 is the second invariant of the deviatoric 
stress tensor and kc are Drucker Prager constants, which are related to the Coulomb’s 
material constants (coefficient of friction /i and cohesion Iq) (see e.g. Bui et ah 2008). For 


the implementation in our SPH code, we use ^/J 2 as a measure of the stress state, and we 


define a pressure dependent yield strength F) for intact rock following Collins et ah (2004): 


Y, = Y. 


RiP 


1 + RiP / {Ym — Yq) 


( 3 ) 


where Yq is the shear strength at P = 0 and Ym is the shear strength at P = 00 and Ri the 
coefficient of internal friction. As in the previous model, F) is temperature depended: 


Y^Y 1 


u 




( 4 ) 


where u is the specific internal energy and Umeit the specific melting energy. 


In the original model by Benz and Asphaug (1994, 1995), fully damaged material was 
treated as strength-less (’’fluid”). As we shall see in section 2.3.1, this simplification leads 


to reasonably accurate results in the case of disruptive collisions between initially intact 
bodies. However, in may situations it is important to take into account the friction in the 
modeling of fully damaged (i.e, granular) material. This is certainly the case when we want 
to study collisions between rubble pile like, granular bodies, or to study the finally shape of 
a body after an impact, and it is also expected to be important in the cratering regime of 
impacts. To model fully damaged rock (damage D = 1), which includes granular material 
in general, we use a yield strength 

Yd = RdP (5) 


where Rd is the coefficient of friction of the damaged material (Collins et ah, 2004). Note 


that Yd is limited to Yd <Yi. In case the modeling starts with a intact or partially damaged 
material, a smooth transition between the criterions (|^ and ([^ is used: 


Y = {D- 1)R + DYd 


( 6 ) 


where Y is limited to H < Hj. 
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If the measure of the stress state \/J 2 exceeds Y, the components of the deviatoric stress 
tensor are reduced by a factor Y/\fj 2 - The yield strength given by replaces the von 
Mises yield criterion used in the previous SPH models. 

Finally, to compute the accumulation of damage D, a tensile brittle failure model based 
on the one of Grady and Kipp (1980) is used 


see 


Benz and Asphaug 

(1994, 

1995 

); 

Jutzi 


et ah (2008). Damage accumulates when the local tensile strain e* reached the activation 
threshold of a flaw. Note that e* is obtained from the maximum tensile stress (j\ after a 
principal axis transformation. 

The friction model described above (equation]^ assumes a constant friction coefficient. 
However, Jop et al. ( |2006 ) use a model with a rate dependent friction coefficient to reproduce 
experiments of dense granular flow. Their 3D model is based on a Drucker-Prager like yield 
criterion and a friction coefficient which is a function of the inertial number, defined as 


/ = miiPip, 


\ 0.5 


(7) 


where P is the isotropic pressure, I 7 I the second invariant of the strain rate tensor 7 jj, d is 
the particle size and ps is the particle density. Using I, the following law for a strain rate 
dependent friction coefficient was proposed 


p(/) — Ps Y {p2 — Ys)/{Io /7 + 1) 


( 8 ) 


where ps is the critical value of the friction coefficient at zero shear rate, p 2 is the limiting 
value at high J, and Jq is a constant. We implemented the relation ([^ in the SPH code by 
replacing the constant coefficient p^ in equation (|^ by p{I). However, as we shall see in 
section 2.3.1, it is not clear whether or not a strain rate depend friction coefficient allows to 
better reproduce the behavior of dense granular material. For this reason, our model uses a 
constant value 

+(7) = Pd (9) 

unless indicated otherwise. 


2.3. Tests 

In this section, we present two test cases where we compare our model with experimental 
results. In the hrst case, we test the ability of our numerical tool to deal with fully fragmented 
(i.e., granular) material. For this, the collapse of a cliff of granular material is simulated. In 
the second case, we test the updated strength models by a comparison to laboratory impact 


experiments using porous gypsum targets (Okamoto and Arakawa, 2009) 


2.3.1. Cliff-collapse 

We model the collapse of a cliff of granular material and we use the results of experiments 


of Lajeunesse et al. (2005) to verify our friction model. Note that Holsapple (2013) performed 


a detailed study of the same problem using a continuum CTH code and Mohr-Coulomb 
and/or Drucker-Prager models. In the cliff collapse problem (see Holsapple 2013 for a 


detailed description), a granular material is initially constrained in a rectangular region 
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with height Hq and length Lq. The width is assnmed to be large and assnmed to be of no 
conseqnence. At the time t = 0, the wall on one side is removed and the material begins 
to flow down dne to downward gravity g. Eventnally, the flow ceases in a hnal rnn-ont 
conhgnration with a with maximnm height H and length L which depend on the properties 


(angle of friction) of the grannlar material. As it was fonnd experimentally (e.g. Lajennesse 
et ah 2005| ) as well as by a scaling analyses (Holsapple, 2013), the hnal scaled prohle only 
depends on the initial height to length ratio and the angle of friction. The gravity and 
the height of the cliff only affect the time and length scales of the problem. This allows to 
model the problem nsing artihcially large dimensions to avoid problems related to time-step 
restrictions (see Holsapple 2013). In the 3D calcnlation presented below, the dimensions 
(inclnding the particle size d) are increased by a factor ~ 10®. 

Fignrej^shows a SPH code calcnlation of a cliff collapse compared with the resnlts of the 
experiments by Lajennesse et ah (2005). In the experiments, glass beads of different sizes. 


and various initial height to length ratios were used. When scaled by the characteristic length 
To and time r = Ho/g, it was found that the prohle curves of the different experiments 
are barley distinguishable. These measured curves are represented by the blue line in Figure 
at f = r = ^yil(jg and f —)■ oo (hnal prohles). 

The results of our simulation using two diherent laws for /i(/) (Equationsandare also 
shown in Figure [Ij The parameters used in the hrst case are fig = tan(20.9), fi 2 = tan(32.76) 


and Iq = 0.279 ( |Jop et al.f 2006), and we use = /is in the second case (constant friction 
coefficient). As it can be seen, both models reproduce very well the experimental results 
and there is no large diherence between the model with a constant /i^ and a rate dependent 
/i(/). The reason is that in the clih collapse problem studied here, the inertial number stays 
small and therefore /i(/) ~ /is = /id- Our results suggest that in the how regime investigated 
by this kind of experiment, both models work equally well, which also means that the global 
outcome of such events is well reproduced by using a single parameter /i^. 

2.3.2. Impact experiment 

In Jutzi et ar| (2009), the strength and fracture model by Benz and Asphaug (1994, 1995) 
in combination with a porosity model (Jutzi et ah, 2008) was compared to laboratory impact 


experiments involving porous pumice. Here, we present a similar comparison test of our 


updated strength model using the results of a study by Okamoto and Arakawa (2009), who 


conducted impact experiments using porous gypsum spheres and a wide range of specihc 
energies (from 3x10^ J/kg to 5x10^ J/kg) and investigated the resulting fragment mass 
distributions. We use run 1 and 2 of this study with the initial conditions indicated in Table 

1 . 

As far as possible, we use material parameters which were measured for the target ma¬ 
terial. For the crush-curve, we use a quadratic form (see Jutzi et ah||2008 ) with Pg = 1x10® 


dyn/cm^ and Pg = 4x10® dyn/cm^, which gives a reasonable ht of the measured crush-curve 
of the gypsum material used in the experiments (Okamoto and Arakawa, private commu¬ 
nication). Further more, the initial slope of the crush-curve in the elastic regime is chosen 
to match the measured longitudinal speed of sound of ~ 2 km/s (Okamoto and Arakawa 


2009). The Weibull parameters were varied in the simulations in order to obtain a match of 


6 



































the largest remnant for one experiment (run 1). For the best £t modeQ m = 9.5 and k = 
8x10^^. For the other runs, the same values of m and k were used. Note that the minimum 
strain threshold for failure resulting from these parameters is about 3x10“® which, multi¬ 
plied by the Young modulus of 5.3x10^^ leads to an equivalent tensile strength of ~ 1.5 
MPa. This value is slightly smaller, but of the same order as the measured tensile strength 
of ~ 2.5 MPa (Okamoto and Arakawa, 2009). In all cases, we use a coefficient of friction 


(damaged material) = 0.55, a coefficient of internal friction (intact material) of /ij = 1.5 
and a limiting yield strength of Ym = 3.5 GP. 

The results of our simulations compared to the experiments are presented in Figure 
1^ The experimental mass distribution is well reproduced in both cases. In Figure the 
results of two simulations of run 1 using different yield criterions are compared. Our results 
indicate that, as long as the same limiting yield strength is used, the Drucker-Prager 
and the originally used von Mises yield criterion lead to very similar results. Besides the 
fragment size distribution, also the antipodal velocities of the fragments was measured. In 
the simulations, we obtain fragment velocities at the antipodal side varying between 4-5 
m/s, which is in good agreement with the experimental values of ~ 5 m/s for run 1 and 2. 
At the center of the antipodal side, we obtain in the simulation also some fully damaged 
material which is ejected with slightly higher speeds (~ 10 m/s). 


3. Catastrophic disruptions 

The catastrophic disruption energy threshold was studied in the past for targets with 
various properties and structures using different previous versions of SPH codes (see Section 
[^. Here, we present results of a new study of using the updated material models 


(see Section 2.2) and investigating the relative effects of friction, porosity and strength. 
Self-gravity is included throughout these simulations using a grid-based gravity solver. To 
identify the largest fragment formed by reaccumulation of smaller pieces, we use an iterative 
procedure based on energy balance (Benz and Asphaug 1999). While this method does 
not allow to compute the whole size distribution of fragments, it was shown to be accurate 


in determining the largest reaccumulated fragment (Jutzi et al. 
signihcant size {Mir/Mfot > 10 - 20%). 

3.1. Comparison to previous results 


2010), as long as it has a 


As a first step, we want to reproduce the results by Benz and Asphaug (1999) for solid 
basalt targets using the updated strength model. We find that, when starting with a intact 
homogenous solid target, both the old (von Mises) and new (Drucker-Prager like) strength 
models lead to the same (within a few %) results in terms of the catastrophic disruption 
threshold as long as the same limiting yield strength Ym is used. Note, however, that 
this is true only in the case where we start with a intact solid target (as it was done in the 


Benz and Asphaug 1999 study). 


^The values for m and k found in this study are different from the ones obtained by Benz and Asphaug 


(1994 19951, because a different kind of material is considered here. 
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3.2. Investigation of various target properties 

Our newly implemented material models allow us to investigate targets with various 
properties. In particular, thanks to the friction model, we can now also study collisions 
among initially fully or partly damaged targets (e.g, granular or rubble pile bodies). In order 
to systematically study the relative effects of friction, porosity and cohesion, we investigate 
the disruption threshold Qf, as a function of radius for spherical bodies with the following 
properties: 


1. no friction, no crushing, no cohesion (purely hydrodynamic) 

2. friction included, no crushing, no cohesion 

3. friction and crushing included, no cohesion 

4. friction, crushing and cohesion included (an intact porous body, similar to pumice) 

In all cases, we start with the same initial density p = 1.3 g/cm^. Targets labeled ’no 
crushing’ also have the same initial density, but there is no crushing of pore space (i.e., 
the initial distention is set to Oq = 1). In the cases that include crushing, a crush-curve 


with parameters for pumice is used (|Jutzi et ah 
compacted 


2009|) which means that the bodies are 
In all the cases, we use a 


(i.e., pore space is decreased) during the collision, 
coefficient of friction (damaged material) pd = 0.8 and a limiting yield strength of Ym = 3.5 
GPa. In the case with cohesion (target 4), a coefficient of internal friction (intact material) 
of Pi = 1.5, and a cohesion Yq = 100 MPa are used. 

For the simulations presented here, we use an impact velocity of 3 km/s and an impact 
angle of 45°. In Figure the outcome of a collision using a target with radius of Rt = 
100 km, and an impactor with Rp= 27 km is shown at the time t = 800s for the cases 
1-4. As it can be seen in the cross sections shown in Figure the degree of disruption and 
consequently the mass of the largest remnant Mir after the collision strongly depend on the 
target properties. The ’hydrodynamic’ target experiences the highest degree of disruption 
and the resulting largest remnant is the smallest {Mir/Mtot = 0.10). In the case 2 where 
friction is included, the degree of disruption is much smaller and Mir/Mtot = 0.67. In the 
cases 3 (with friction and crushing) and 4 (with friction, crushing and cohesion), the target 
is signihcantly compacted by the collision and the density is increased accordingly. In both 
cases, Mir/Mtot = 0.82, which is larger than for the other targets. The outcome for runs 
3 and 4 are very similar, which indicates that for disruptive collisions at this size scale, 
the effects of cohesion and tensile strength are negligible. On the other hand, porosity 
(energy dissipation by compaction) and friction (when starting with a damaged body) are 
very important. Figure shows the disruption threshold Q/id for the for targets 1 - 4 as a 
function of the combined target radius Rci (defined in section [^. The curves confirm 
the results discussed above for the Rt = 100 km case. The hydrodynamic target has a 
by far the lowest QfiQ. The disruption threshold is significantly (5-10 times) higher when 
friction is included and is further increased (by a factor of 2-3) when the energy dissipation 
by compaction (pore crushing) is taken into account. The curves for the targets 3 and 4 are 
almost exactly the same for target sizes Rt>l km. For smaller targets, cohesion and tensile 
strength start to affect the outcome and at even smaller scales, the largest fragments in the 
case of the cohesive targets are intact fragments rather then purely gravitational aggregates, 
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which explains the deviation of the cnrves of the targets 3 and 4 at small sizes (strength 
regime). As it was fonnd in previous studies, the increases with decreasing size for 
cohesive targets in the strength regime due to the size dependence of the tensile strength. 
In the Figure we also plot one point representing the Q*pid for the Rt = 50 km rubble 
pile bodies used in the study by Benavidez et ah (2012). In this study, the rubble pile 


targets were constructed by hlling the interior of a 100-km-diameter spherical shell with an 
uneven distribution of solid basalt spheres having diameters between 8 km and 20 km. It 
is important to note that this study was performed using an old version of the SPH code 
which does not include friction and which uses the von Mises yield criterion. As shown in 
Figure the rubble pile targets used in the Benavidez et ah (2012) study have a very low 


Q*jid which even falls below the curve of our hydrodynamic targets. 


4. Discussion and ontlook 

In this paper, we presented recent improvements of the SPH technique concerning the 
modeling of the disruption of strength dominated bodies. The updated models are able to 
reproduce the results of laboratory experiments very well. As for the modeling of catas¬ 
trophic disruptions of asteroids, a comparison between old and new strength models shows 
no signihcant deviation in the case of targets which are initially non-porous, fully intact and 


have a homogenous structure (such as the targets used in the Benz and Asphaug 1999 study). 
In this case, the crucial parameter is the limiting yield strength Ym (see also Leinhardt and 


Stewart 2009), and the details of the strength model (e.g., pressure dependent vs. pressure 
independent yield strength) do not play an important role. At large scales (> 1 km) the 
bodies get fully damaged in the disruptive collisions investigate here. It was found that in 
this case, the effect of friction in the post-impact flow (after the shock wave has damaged the 
material) is negligible. However, it is important to point out that this is not true in many 
other cases (e.g., initially partly or fully damaged targets, rubble-pile structures, etc.) in 
which it is crucial that the damaged material still has a pressure dependent shear strength 
(i.e, there is friction). As our investigations of the disruption threshold for different target 
properties show, an initially (pre-impact) fully damaged target modeled without friction or 
compaction has a Q*jid which is signihcantly (5-10 times) smaller than the for the same 
target but where friction is taken into account. Interestingly, the rubble-pile targets used in 


the study by Benavidez et ah (2012) (who used an older version of the SPH code) seem to 


have a Q*^ which is very low and even lies below our Q*^ curve for the purely hydrodynamic 
bodies. Our results therefore indicate that these bodies represent something like a porous 
fluid rather than real rubble-pile bodies, which we believe should have a Q*jj lying between 
the curves of the cases 2 and 3 in Figure As our investigations show, when the effect of 
the energy dissipation due to compaction (pore crushing) is taken into account, the targets 
become stronger is increased by a factor of 2-3)R On the other hand, cohesion is only 


^Note that Jutzi et al. (20101, in a comparison between porous pumice and solid basalt targets, found 
that at large scales (gravity regime), the non-porous targets have a higher This different behavior is 

due to differences in density and shear strength (1m) 
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important at small scales (< 1km) and has an negligible effect at larger scales. Our results 
therefore indicate that in the gravity regime, the crucial parameters are the crushing prop¬ 
erties (crush-curve), friction (when starting with partly or fully damaged bodies) and the 


limiting yield strength Ym (the effect of varying Ym was already investigated by Leinhardt 


and Stewart 2009). Cohesion and tensile strength do not play an important role in this 


regime. 

Our results clearly confirm that collisions between bodies of a few 100 km diameter can 
not be treated by just using ’hydrodynamic objects’, even if they are fully damaged without 
cohesion. However, at a certain size, bodies are expected to be fully gravity dominated (due 
to the increase of overburden pressure with size); the transition to this regime will be subject 
of a future study. 

The results presented here, as well as those from an ongoing, more general study (where 
a much larger parameter space will be covered) will be used to improve existing scaling laws 
for the outcome of collisions (e.g. Leinhardt and Stewart 2012). 
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Figure 1: Collapse of a granular cliff. The height (y-axis) is normalized by the initial height Hq and the 
length (x-axis) is normalized by the initial length Lq. The experimental results by Lajeunesse et al. (20051 
(light blue) are compared to the results of SPH code calculations using the model for dense granular flow 
by Jop et al. ( |2006 ) (solid red line) and a Drucker-Prager friction model with a constant friction coefficient 
(dotted blue line). The profiles are shown at the times t = t — ^H^jg and t —>■ oo (final profiles). 
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Figure 2: Cumulative fragment size 
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simulation (red). The targets consist of porous gypsum; the impact 
and rnn 2 (bottom) of the experiments described in Okamoto and 
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Figure 3: Same as Figurej^ but only run 1 is shown. In the SPH simulations, the Drucker-Prager (DP) and 
von Mises (VM) strength models are compared. 
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Figure 4: Cross-section of SPH simulations of collisions between a target with a radius of Rt = fOOkm and 
projectile of Rp = 27 km with a relative velocity of 3 km/s and a 45° impact angle. Four different targets 
f - 4 (as indicated at the top) are investigated. The outcomes in terms of the degree of disruption, size of 
the largest remnant {Mir/Mtot) and density increase strongly depend on the target properties. 
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Figure 5: Catastrophic disruption threshold as a function of the combined radius Rci for various target 
properties. The curves are the result of new SPH code calculations. The black point corresponds to 
the value of estimated from simulations of ”rubble-pile” collisions by Benavidez et al. (20121. 
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Table 1: Initial conditions and results of impact experiments ( |Okamoto and Arakawa 2009). mi/Mt is the 
largest fragment mass normalized by the original target mass (Mt). Vi, Q and 14 are the impact velocity, 
the specific energy and the antipodal velocity, respectively. 


Run number Mt (g) Target diameter (mm) mi/Mt Vi{km/s) Q (J/kg) Vaijn/s) 

No. 1 TTI 28)8 0.377 3)33 3.21 x 10^ 5.08 

No. 2 8.44 25.5 0.281 3.30 4.51 x 10^ 5.00 
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